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Li diffusion in Lia-CoC^ probed by muon-spin spectroscopy 
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The diffusion coefficient of Li + ions (Du) in the battery material Li^CoC^ has been investigated 
by muon-spin relaxation (/x + SR). Based on the experiments in zero- field and weak longitudinal- fields 
at temperatures up to 400 K, we determined the fluctuation rate (u) of the fields on the muons due 
to their interaction with the nuclear moments. Combined with susceptibility data and electrostatic 
potential calculations, clear Li + ion diffusion was detected above ~ 150 K. The Dm estimated from 
v was in very good agreement with predictions from first-principles calculations, and we present the 
/i + SR technique as an optimal probe to detect Du for materials containing magnetic ions. 

PACS numbers: 76.75.-M, 66.30.H-, 82.47.Aa, 82.56.Lz 
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In spite of a long research history on lithium insertion 
materials for Li-ion batteries [l[, e.g., LiCoC>2, LiNiC>2, 
and LiMn2C>4, one of their most important intrinsic phys- 
ical properties, the Li + ions diffusion coefficient (Djj), 
has not yet been determined with any reliability. Al- 
though Li-NMR is, in general, a powerful technique to 
measure Z?Li for non-magnetic materials, it is particularly 
difficult to evaluate Z?Li for materials containing mag- 
netic ions, because the magnetic ions induce additional 
pathways for the spin-lattice relaxation rate (1/Ti), re- 
sulting in huge 1/Ti compared with that expected for 
only the diffusive motion of Li ions. 

Such difficulty was clearly evident in the \/T\(T) curve 
for LiCoC>2 and LiNiC>2 @, and, for that reason, Th^ 
was instead estimated from the Li-NMR line width [J]. 
However, since the line width, i.e., the spin-spin relax- 
ation rate (I/T2) is also affected by the magnetic ions, the 
D Li obtained by Li-NMR for LiCo0 2 (= 1 x KT 14 cm^" 1 
at 400 K) is approximately four orders of magnitude 
smaller than predicted by first-principles calculations [f| . 
Since lithium insertion materials always include transi- 
tion metal ions, in order to maintain charge neutrality 
during the extraction and/or insertion of Li + ions, it 
is consequently very difficult to determine Z?Li for these 
compounds unambiguously by Li-NMR. 

On the other hand, the chemical diffusion coefficient 
(D^ cm ), which is measured under a potential gradient, 
has also been determined by electrochemical measure- 
ments. Note that the relationship between Du and 
£,chem jg giycn by £>chcm = Q£) U: where 9 is a ther- 
modynamic factor. The magnitude of D^ cnl is, however, 
known to be very sensitive to the measurement system, 
e.g., the electrolyte as well as the compositions of the pos- 
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itive and negative electrodes. As a result, the reported 
^chem for Li^CoOz ranges from 4 x 10~ 8 to 10~ 10 enr^s" 1 
for powder samples [1, M H, H, GU an d fr° m 2-5 x 10~ n 
to 2 x 10" 13 cn^s" 1 for thin films [U GJ, GJ, at am- 
bient T. In order to profoundly understand the physics 
behind the operation principle of battery materials, it is 
imperative to have a reliable probe to measure Du for 
all the components of the battery as a function of both 
Li content as well as T. This is at present a key issue for 
further development of Li-ion batteries, and in particular 
for future fabrication of solid-state batteries. 

In contrast to NMR, the effect of localized moments in 
a paramagnetic (PM) state on the // + SR signal is very 
limited at high T, because the PM fluctuation is usually 
too fast to be visible by /i + SR. Furthermore, although 
/i + SR is very sensitive to the local magnetic environment, 
whether it is due to electronic or nuclear spins, an elec- 
tronic contribution is, in principle, distinguishable from 
a nuclear contribution by weak longitudinal field mea- 
surements. In addition, since the /i + SR timescale is dif- 
ferent from the NMR's one, ^ + SR is expected to provide 
unique information on nuclear magnetic fields. Indeed, 
/i + SR studies of LIeCoC^ at low T indicate that /i + s 
feel a nuclear magnetic field caused by Li and 59 Co even 
at 1.8 K 0, [l6|]. This offers a possibility to determine 
TJ L i of Lij;Co02 by /i + SR, if /i + s locate in the vicinity 
of the 2 ~ ion and make a stable /j, + -0 2 ~ bond in the 
LizCoC^ lattice. Here, we report our initial work on 
Lia;Co02 (x=0.73 and 0.53) to estimate their Du and 
establish yU + SR as a novel tool to probe Li-ion diffusion. 

A powder sample of LiCo02 was prepared at Osaka 
City University by a solid-state reaction technique us- 
ing reagent grade LiOH-H20 and C0CO3 powders as 
starting materials. A mixture of the two powders was 
heated at 900°C for 12 h in air. Powder X-ray diffrac- 
tion (XRD) analysis showed that the LiCo02 sample was 
single phase with a rhombohedral system of space group 
R3m (oh = 0.2814 nm and ch = 1.4049 nm in hexag- 
onal setting). The Li-dcficicnt samples were prepared 
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FIG. 1: (Color online) ZF- and two LF-/i + SR spectra for 
Lio.73Co02 measured at (a) 100 K and (b) 225 K. The mag- 
nitude of LF was 5 and 10 Oe. Solid lines represent the fit 
result using Eq. |T}. 
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by an electrochemical reaction using Li|LiPFg-ethylene 
carbonatc-dicthyl carbonate|LiCo02 cells. The LiCo02 
powder was pressed into a disc with 15 mm diameter 
and 0.4 mm thickness, and the disc was then used as a 
positive electrode. The Li x Co02 disk was removed from 
the cell in a glove-box and packed into a sealed powder 
cell just before the /x + SR measurement. Their structures 
were subsequently confirmed by powder XRD, and, fi- 
nally, their compositions were checked by an inductively 
coupled plasma atomic emission spectral analysis. The 
above procedure is essentially the same as that of our 
previous /x + SR work on LLjCoG^ [HI and LLjNiO^ [l7| . 

The /i + SR spectra were measured at the ARGUS sur- 
face muon beam line of the RIKEN-RAL Muon Facility 
at ISIS in the UK using a liquid-He flow type cryostat 
in the T range between 10 and 400 K. The experimental 
techniques were described elsewhere [Til ], x was mea- 
sured using a SQUID magnetometer (MPMS, Quantum 
Design) in the T range between 5 and 200 K under a 
magnetic field of H = 100 Oc. 

Figure Q] shows the zero field (ZF-) and longitudinal 
field (LF-) yU + SR spectrum for the Lio.73Co02 sample 
obtained at 100 and 225 K. At 100 K, the ZF-spectrum 
exhibits a typical Kubo-Toyabe (KT) behavior with a 
minimum at t ~ 6 /xs, meaning that the implanted muons 
see the internal magnetic field (-ffint) due to the nuclear 
magnetic moments of 7 Li, 6 Li and 59 Co. The applied LF 
clearly reduces the relaxation rate, i.e., the time slope, by 



FIG. 2: (Color online) For Lio.73Co02 and Lio.53Co02 re- 
spectively, we show T dependences of (a, b) field distribution 
width (A), (c, d) field fluctuation rate (y), (e, f) susceptibil- 
ity (x), and (g, h) the relationship between \og{v) and 1/T. 
A and v were obtained by fitting both ZF- and LF-spectra 
using Eq. ([T]). \ was measured in both field cooling {FC) and 
zero field cooling (ZFC) mode with H=W0 Oe. In (e, f), 
the x data [l(J measured in FC mode with H — 10 kOe for 
Lio.75Co02 and Lio.52Co02 were also plotted for comparison. 
The straight lines in (g) and (h) show the activated diffusive 
behaviour discussed in the text. 

decoupling H mt . Although the ZF-spectrum still shows 
KT behavior at 225 K, the relaxation rate is smaller than 
at 100 K. 

In order to estimate the KT parameters precisely, 
the ZF- and two LF-spectra were fitted simultaneously 
by a combination of a dynamic Gaussian KT function 
[G KT (A, v, t, Hlf)] and an offset signal from the frac- 
tion of muons stopped mainly in the sample holder, which 
is made of high-purity aluminum; 

A P LF (t) = A KT G DGKT (A,v,t,H LF ) + A BG (1) 

where Aq is the empirical maximum muon decay asym- 
metry, Akt and Abg are the asymmetries associated 
with the two signals. A is the static width of the lo- 
cal field distribution at the disordered sites, and v is 
the field fluctuation rate. When v = and i?LF = 0, 
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FIG. 3: (Color online) The variation of electrostatic potential 
(4>e) along the a-axis at b = b m and c = (on the Li plane) 
and 6 = and c = 0.022 (1 A away from the O 2- ions) in 
the hexagonal lattice of Li^CoO^ with x — 0, 0.25, 0.5, 0.75, 
and 1 (from top to bottom). Here, b m denotes the 6, at which 
electrostatic potential exhibits a minimum. The right panels 
show the corresponding distribution of (/>e in the Li plane. 
Arrows on the 6-axis in the right panels represents b m . 

G DGKT (t, A, v,Hi, F ) is the static Gaussian KT function 
GfJ(t, A) in ZF. At first, wc fitted all the ZF-spectra 
using common Akt and Abg in Eq. |T|). The "global 
fit" provided that A KT = 0.164056±0.000011 (0.16889± 
0.00018) and ^ B g = 0.06350 ± 0.00002 (0.0692 ± 0.0002) 
for Lio.73Co02 (Lio.53Co02). Then, using the obtained 
Akt and Abg > one ZF- and two LF-spectra were global- 
fitted using common A and v at each T. 

Figure [5] shows the T dependencies of both A and v for 
the two samples together with x measured in a 100 Oe 
magnetic field. For L10.73C0O2, A is almost constant for 
5 < T < 250 K, indicating that the are most probably 
stable in the crystal lattice until ~ 300 K. A then rapidly 
decreases, but levels off again for T > 325 K. The v(T) 
curve is almost T-indcpcndent up to 150 K, starts to 



increase at ~ 150 K, and exhibits a maximum at 275 K. 
Above 275 K, v decreases to 0.7 x 10 6 s" 1 at ~ 325 K, and 
finally becomes almost T-independent above 350 K. The 
increase in v between 150 and 275 K is well explained by 
a thermal activation process [Figs.[2]Jg) and[2jh)], which 
signals the onset of diffusive motion of either Li + or /i + 
above 150 K. The clear decrease in A at ~ 300 K also 
suggests an additional diffusion of Li + or /j + . 

The x(T) curve exhibits a small anomaly around 150 K 
with a thermal hysteresis of ~ 10 K, while there is no 
indication of any magnetic anomalies in the T range be- 
tween 200 and 300 K [Figs.^c) and Hp)] [H,[I|. This 
suggests that the change in the ^ + SR parameters around 
150 K is caused by an intrinsic change in Li^CdO^, but 
the change around 300 K is visible only by /i + SR. The in- 
crease in v above 150 K is, thus, most unlikely due to [x + 
diffusion but it is rather due to Li + diffusion, i.e., either 
a freezing of the Li + motion or an order-disorder tran- 
sition of the Li + ions occurs below around 150 K. This 
is also supported by a recent 7 Li-NMR experiment [lj|, 
in which the NMR line width -vs.-T curve for Lio.6Co02 
exhibits a step-like decrease with T around 150 K by 
motional narrowing due to Li + diffusion. Since such dif- 
fusion naturally increases a local structural symmetry, 
it is reasonable that A slightly decreases with T around 
150 K. On the other hand, both Li + and ^ + are inferred 
to be diffusing above 300 K, resulting in the large de- 
crease in A caused by motional narrowing. Actually, be- 
cause A < OAv above 300 K, Eq. ([T]) is roughly equivalent 
to an exponential relaxation function [exp(— At)], and it 
is difficult to estimate A and v precisely at high T. 

The result for L10.53C0O2 sample is very similar 
to that of Li .73CoO2, although AT->o(Lio.53Co02) < 
At^o(L1o.73Co02) due to the decrease in the number 
density of Li + ions, as reported previously [HI, Also, 
the magnitude of v of Lio.53Co02 is smaller than v of 
Lio.73Co02 in the whole T range measured, but the v(T) 
curve for both samples show a clear increase with T above 
150 K and a maximum around 300 K. 

In order to predict the muon site(s) and to confirm 
the reliability of the above assumption that Li + ions dif- 
fuse above 150 K whereas /i + diffuse only above 300 K, 
wc performed electrostatic potential (0e) calculations for 
the LizCoO^ lattice using a point-charge model and the 
program DipElec [201 ]. As seen in Fig. [31 the site in 
the vicinity of the 2 ~ ions is more stable for fj, + than 
the site in the Li plane for the whole x range between 
1 and 0. This means that /i + 's are bound to the O 2- 
ions so as to make a stable fi + -0 2 ~ bond in Li^CdO^. 
This is a common situation in oxides, as for example in 
the case for the high-T c cuprates [2l|. Since ^t + s are 
assigned as an ideal point charge, such /i + -0 2 ~ bond 
should be purely ionic. In fact, dipole field calcula- 
tions for the site in the vicinity of the 2 ~ ions pro- 
vide that A calc = 0.43 x 10 6 s^ 1 (0.35 x 10 6 s" 1 ) for 
Lia;Co02 with x = 3/4 (1/2). Furthermore, A ca i c is 
found to be comparable to A measured at low T in the 
whole x range for Li^CoO^ [3, if wc consider the re- 
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FIG. 4: (Color online) (a) Possible jump paths for Li ions. 
Broken arrows represent the direct jump to nearest (vacant) 
Li site (path No. 1), whereas solid arrows the jump to an 
interstitial site in the center of the oxygen tetrahedron (path 
No. 2). (b) The relationship between 7?Li and x in Li^CoO^ 
as extracted from our /i + SR experiment (open circles). Solid 
and dashed lines represent the predictions by first-principles 
calculations at T = 300 and 400 K, respectively, when 
the effective vibration frequency is a typical value (10 13 s _1 ). 
Sharp minima in the predicted curve (at x — 1/3 and 1/2) 
are caused by Li-ordering. The 7 Li-NMR result 0, Qj| for 
T = 400 K is also plotted (solid dot) for comparison. 



duction of A by the electric field gradient effect on the 
nuclear moments with I > 1 [22l , l23j . This suggests that 
the point-charge model is acceptable for determining the 
muon site(s) in Li^CoC^. As a result, it is clarified that, 
as T increases from 5 K, the Li + ions start to diffuse 
above 150 K (= T^ 1 ) and then the fi + diffuse above 
300 K (= Tj 1 ), in spite of the mass difference between 
fi + and Li + (m L j+ /m fl + ~ 63) because the muons form 
a hydrogen-like bond with oxygen. 

Finally, we estimate Z?Li using the obtained fluctuation 
rate v as directly measuring the jump rate. Figure Ufa) 
shows the possible jump paths for the Li ions. That 



is, the direct jump to the nearest (vacant) Li site (path 
No. 1) and the jump to the interstitial site in the center 
of the oxygen tetrahedron (path No. 2). Assuming that 
v corresponds to the jump rate of the Li ions between 
the neighboring sites, Du is given by [24| ; 
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where TV* is the number of Li sites in the i-th path, Z v ^ is 
the vacancy fraction, and Sj is the jump distance. Here, 
we naturally restrict the path to lie in the c-plane, i.e., 
along the 2D channel, because it is most unlikely that 
the Li ions jump across the C0O2 plane to an adjacent 
Li plane. Therefore, N± = 6, N2 = 3, Sx is equivalent to 
the a-axis length, s 2 = a/y/3, Z Vy i = 0.27 for L10.73C0O2 
(0.47 for L10.53C0O2), and Z v> 2 = 1- As a result, we ob- 
tain D Li = (7±2) x 10~ 10 cm 2 /s [(2.5±0.8) x 10~ 10 cm 2 /s] 
for L10.73C0O2 [Lio. 53 Co0 2 ] at 300 K. Here, ^(300 K) for 
L10.73C0O2 was estimated from the extrapolation of the 
linear relationship between log[f] and T -1 [see Fig.[2jg)]. 
The estimated Du is found to be very consistent with 
the prediction by first-principle calculations [f|, as seen 
in Fig. [31(b) . Note that the jump paths used in Eq. @ 
are the same to those for the first-principle calculations. 
This means that there is no ambiguous factor for estimat- 
ing Du by yU + SR. Since /i + SR detects v ranging from 
~ 0.01A to ~ 10A, it is applicable for materials with 
D u = 10~ 12 - 10~ 9 cm 2 /s, when TV = 10, Z v = 1, 
s = 1 nm, and A = 0.1 x 10 6 s^ 1 . 

In conclusion, we have been able to determine the Li 
diffusion coefficient, Du, of Li^Co02 from the fluctuation 
rate of the field experienced by the muons in interaction 
with the nuclear moments of the diffusing ions. The value 
was found to be in good agreement with theoretical pre- 
dictions. Consequently, we would like to suggest /x + SR 
as a novel probe to investigate Li diffusion, especially for 
materials containing transition metal ions. 
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